Intrinsic Spin Hall Effect: Topological Transitions in Two-Dimensional Systems 
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The spin-Hall conductivity in spatially-homogeneous two-dimensional electron systems described 
by the spin-orbit Hamiltonian htt p -a is presented as a sum of the universal part Me/Snfi determined 
by the Berry phase $ = Mn (M is an odd integer, the winding number of the vector fi p ) and a 
non-universal part which vanishes under certain conditions determined by the analytical properties 
of ftp. The analysis reveals a rich and complicated behavior of the spin-Hall conductivity which is 
relevant to both electron and hole states in quantum wells and can be detected in experiments. 

PACS numbers: 73.63.-b, 72.25.-b, 72.25.Pn, 71.70.Ej 



Owing to the spin-orbit interaction (SOI), an electric 
field applied along two-dimensional (2D) electron lay- 
ers can generate transverse spin currents in the absence 
of external magnetic fields. This phenomenon, known 
as the spin Hall effect [1], is at the focus of attention 
in modern physics. The presence of SOI terms in the 
Hamiltonian of free electrons leads to the intrinsic spin- 
Hall conductivity expressed in the universal units e/h in 
the case of weak disorder. The original theoretical pro- 
posal [2] of the universal intrinsic spin Hall effect has 
been based on the Rashba Hamiltonian [3] describing 
the linear in 2D momentum p = (p x ,Py) spin-orbit cou- 
pling due to structural inversion asymmetry. However, 
numerous theoretical calculations [4-8] have proved the 
absence of static intrinsic spin currents for this case. As 
follows from the equation of motion for the spin density 
operator [9], this statement is applicable to any electron 
system described by a p-linear SOI Hamiltonian. The 
situation is quite different in the case of 2D hole sys- 
tems described by the effective p-cubic SOI Hamiltonian 
[10] hp^ — ?ik((7 + p 3 _ + <7_p+), where a± = (a x ± ia y )/2, 
a a are the Pauli matrices, and p± = p x ± ip y . Theo- 
retical studies [11-15] based upon this Hamiltonian have 
confirmed the existence of the intrinsic spin Hall effect. 
The experimentally observed spin Hall effect in 2D hole 
systems [16,17] is likely of the intrinsic origin. 

What makes the systems described by the Hamilto- 

(3) 

nian h p so different from the systems described by p- 
linear SOI Hamiltonians? It is the dependence of SOI on 
the angle tp of the 2D momentum. This dependence is 
characterized by the odd integers known as the winding 
numbers (WN), which are equal to ±3 and ±1 for the 
p-cubic and p-linear Hamiltonians considered above. In 
2D hole systems, owing to the increased WN, the con- 
servation of the spin density is no longer reduced to the 
requirement of zero spin currents, so the intrinsic spin 
Hall effect exists. The role of WN in spin response can 
be also emphasized by considering their influence on the 
collision-mediated spin-charge coupling term known from 
the Kubo formalism as the vertex correction [4]. If the 



scattering is symmetric (caused by the short-range po- 
tential), the vertex correction for the WN ±3 is zero, 
since it is given by the angular average of the product of 
the charge current operator by the SOI Hamiltonian. In 
contrast, for p-linear SOI the vertex correction is always 
nonzero and leads to nonexistence of spin currents. 

It is important to realize that the consideration of SOI 
Hamiltonians containing the terms with WN either ±1 
or ±3 is not sufficient for description of the spin response 
in 2D systems. The coexistence of SOI terms with WN 
±1 and ±3 in semiconductor quantum wells is rather a 
rule than an exception. For example, this is the case of 
conduction-band electrons in the quantum wells made of 
noncentrosymmetric semiconductors [18] at high electron 
densities, when both p-linear and p-cubic Dresselhaus 
terms are important. The aim of this Letter is to find 
out the general properties of the intrinsic spin currents 
for the systems described by the SOI Hamiltonians con- 
taining an arbitrary mixture of terms with different WN 
and to establish relevance of such a consideration to both 
electron and hole states in quantum wells. 

The starting point is the free-electron Hamiltonian in 
the momentum representation: 



H p = e p + h p , h r 



(i) 



where s p is the kinetic energy (isotropic but not necessar- 
ily parabolic) and Jl p is an arbitrary vector antisymmet- 
ric in momentum. The 2x2 matrix SOI term h p describes 
both 2D electrons and 2D holes (since the 4-fold degener- 
acy of the Tg, valence band is lifted in quantum wells, the 
2D holes are quasiparticles with two spin states). The 
calculations are based on the quantum kinetic equation 
for the Wigner distribution function [19] which is a 2 x 2 
matrix over the spin indices. Searching for the linear re- 
sponse to the applied electric field E in the stationary 
and spatially homogeneous case, one can write the dis- 
tribution function in the form fp 9 * 1 + / p , where f p is 
the non-equilibrium part satisfying the linearized kinetic 
equation 



hp, fp 



eE 



df p eq) 
dp 



= J(f\p). 



(2) 
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The collision integral J describes the elastic scattering, 
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and the spin-orbit corrections [1] to the scattering poten- 
tial are neglected. Considering this integral in the Marko- 
vian approximation and assuming that htt p = h\€l p \ is 
small in comparison to the mean kinetic energy, one can 
expand J in series of f2 p [14,19]. Using the spin-vector 
representation / p = f p + & ■ f p , one gets 



- 2[fi p x f p ] + A p = ^ [*% w \*-v>\Qv> fp). 



(3) 



where w q is the Fourier transform of the correlator of 
the scattering potential, |p'| = |p| is assumed, and tp' is 
the angle of the vector p'. Next, A p is a vector propor- 
tional to E, and m p = \{dp 2 / de p ) is the p-dependent 
effective mass as it enters the expression for the group 
velocity, v p = de p /dp = p/m p . Analytical solution of 
Eq. (3) is possible for short-range scattering potential, 
when W[ p _ p /| ~ w. Then the right-hand side of Eq. (3) is 
written as is p ({ p — f p ), where v p = m p w/h 3 is the scatter- 
ing rate and the line over a function denotes the angular 
averaging. Also, 



he 



p _ ' vp • P--p •■' : p + j E x ^ fe p dp^lp 



de 



v 



p 2 dtp 



, (4) 



where f' E = df e / de p is the derivative of the Fermi dis- 
tribution function f £ , and n is the unit vector normal to 
the quantum well plane. Notice the property A™ = 0. 

Solution of Eq. (3) determines the non-equilibrium 
spin current density q 7 = \ J -^^Tt({&, u 7 (p)}/ p ), 

where u 7 (p) = d(e p + h p )/dp 1 is the group velocity 
in the presence of spin-orbit interaction, { , } denotes 
the symmetrized matrix product, and Tr is the matrix 
trace. The spin conductivity is introduced according to 
q 7 = TiypEp. Based on Eqs. (3) and (4), 



e 



■ J de p f' £p {Tf-[PlxQ^}). (5) 

The vector-functions standing here are defined as angular 
averages: Tf = 2[Pg x (dSlp/dp p)], = P 7 ft P /A p , 

A p = + and Q£ = 2[R p ]~ 1 R p (dn p /d P p), where 
R p is a symmetric matrix with elements R p P = (QpSap — 
Opf2^)/A p . One can find also the induced spin density: 
s = U (W Tr (^/p) = (eh 2 /4-Kw)fds p f' Ep Q$Ep. The 
limit of low temperature [20] is described by the substi- 
tution f = — 5(s p — ef), so the spin conductivity tensor 
is expressed directly through the vector-functions taken 
at the Fermi surface e p = e PF = Ef- 

Equation (5) is valid for arbitrary f2 p . In the quantum 
wells grown along [001] direction in cubic crystals of zinc- 
blende type, the C2v point group symmetry implies 



(6) 



where the polar coordinate representation p = (p, tp) is 
used. Then Tf = (0,0, Tf), V% = {Pp,P^,0), and 



= (Qf ,Qf,0), where T? 



rpyy rp X y _ r TyX 

-Qyy, and Q;y = 



pxx _ pyy pxy _ pyx f)xx 

r p " r p i r p r p i Wp 

— Q v v x ■ The spin currents exist only for z-spins, £ 7/ 3 = 
(0, 0,S 7l g), and there are two independent components 
H xy = —T, yx = Yih and T, xx = —^ yy describing spin-Hall 
and spin-diagonal currents, respectively. The function 
T p y entering E# can be written as 



rpxy _ 
V 



P:<P 



2ttA^ V p - v dp 



(7) 



In the case of zero temperature, using the notations 
T xy = T p y and VL^ = Sl PF ^, it is convenient to write 



mh 



(8) 



where <5£# expresses the contribution of the second term 
in Eq. (5). In the collisionless limit, the formal integra- 
tion in Eq. (7) leads to 

T xy =®, ^ = ^Ida.rg[n+(z)}=n(N -N oo ), (9) 



iO.y, is a function of the complex 



where fl + (z) 

variable z = e ltp , and the contour of integration in the 
complex plane is the circle of unit radius, \z\ = 1. Next, 
iVo and are the numbers of zeros and poles of fl + (z) 
inside this circle (it is assumed that Q + (z) does not have 
branch points). Using the conventional definitions (see 
[21] and references therein) it is easy to identify $ with 
the Berry phase in the momentum space. In the WN 
representation, the function fi + (z) is a polynomial con- 
taining odd powers of z, in the general case, from z~ N 
to z N , assuming that the highest WN involved in i~2 p is 
N. Then = L, where L < N is an odd integer (the 
order of the multiple pole at z = 0), while N takes even 
values from to L + M, where —L<M<N depending 
on the SOI parameters. Therefore, if tt p contains an ar- 
bitrary mixture of terms with different WN up to N, the 
spin-Hall conductivity is 



Me 
8^h 



+ 6Z H , |M| = 1,3,...,JV, (10) 



where M = No — Noo is the acting WN, which describes 
the actual winding of the vector fl p as p goes around 
the Fermi surface, and can be found, in each concrete 
case, from the simple analysis explained above. The cor- 
responding Berry phase is irM. The spin-Hall conduc- 
tivity changes abruptly when the functions Sl x and Qv, 
go through zero simultaneously at certain angles tp. In 
other words, each time when the SOI parameters are ad- 
justed in such a way that the spin splitting 2hfl p at the 
Fermi surface becomes zero at certain p, a topological 
transition occurs: the Berry phase changes by ±2ir. For 
the Hamiltonians with N = 1 including both Rashba 
and Dresselhaus (linear) terms, this effect has been stud- 
ied in the Berry phase approach in Refs. 21-23. In this 
particular case, however, the first term in Eq. (10) is 
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exactly compensated by the second term, and = 0. 
Therefore, the topological transitions essentially require 
the SOI with WN greater than unity. 

The result (10) is exact in the collisionless limit and can 
be viewed as a quantization of the spin-Hall conductivity 
in terms of the WN. In general, this quantization does not 
occur in integer numbers of e/Airh, because <5£# is also a 
discontinuous function of SOI parameters and undergoes 
abrupt changes together with the first term in Eq. (10). 
To show this, it is sufficient to represent P aj as combi- 
nations of the integrals § dz[Vl + (z)]^ 1 , § dz[z 2 ^l + (z)]^ 1 , 
and complex conjugate terms. It is important that such 
a representation allows one to find the general conditions 
for vanishing <5£#: this takes place when either a) all ze- 
ros of fl + (z) are inside the circle H = 1 and No — N x > 3 
or b) the order of the multiple pole at z — is L > 3 
and all zeros of f2 + (z) (if present) are outside the cir- 
cle \z\ = 1. In particular, this means that if the highest 
WN involved in fi p is, in the same time, the acting WN 
(M = N or M = -N at N ^ 1), the spin-Hall con- 
ductivity stays at the universal value Me/8irfi without 
regard to the SOI parameters. If N > 5, can take 
universal values from ±3e/87r?i to ±Ne/8nh. 



7i are the Luttinger parameters in their usual notations, 



leads to the 6„-term with b p = jjLa p . 



Next, the c p - and 



Mp-terms for holes are caused by the bulk inversion asym- 
metry [24]. The Cp-term includes the contribution cthP 
proportional to p [24,25], which should dominate at low 
hole densities. In the general case, especially when the 
structural asymmetry is weak, an adequate description 
of hole states should include all terms in Eq. (11). 
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FIG. 2: Spin-Hall conductivity as a function of density in elec- 
tron (a) and hole (b) systems. The dashed lines correspond 
to the collisionless approximation, v = 0. The solid lines are 
plotted for v = 0.2/3(irh/d) (a) and for v = 0.5|a|/«;| 1/2 (b). 
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FIG. 1: Left: Phase diagram for the SOI of Eq. (11) at 
Up — 0. The regions of fixed Berry phase $ (indicated) 
are separated by the lines of topological transitions (solid). 
Right: Spin-Hall conductivity Tih (solid) and its universal 
part (dash) as functions of b/a at c/a = 1. It is assumed that 
9 In \b p \/d\np = 9 In \c p \/d\np = 1 and d In |o p |/91np =3. 



The simplest case of the SOI with combined WN de- 
scribed by Eq. (11) is realized when c p — u p = 0. One 
finds the analytical expression 



8irh 2b 2 r 2 



[(S-r^-^ + r/r 2 ], (12) 



where r p 



+ M 2 +^ 2 /4][K - fo p ) 2 +^/4]-^ 2 /4, 

i] p = l/2+(l/4)d\n\b p /a p \/dlnp, and all coefficients are 
taken at p = pf- According to the Berry phase analysis, 
T,h = 3e/87r7i at |a| > |6| in the collisionless limit, while 
Eq. (12) gives 



3e 

8^h 



1, a 2 > b 2 

-(l-2 V /3)(a/b) 2 , a 2 < b 2 



(13) 



The most general form of f2 p including WN ±1 and 
±3 for [001]-grown quantum wells is 

ftp = c p sin tp — bp cos ip — u p sin 3^ — a p cos 3ip, 
fl^ — — c p cos ip + bpSinip — u p cos3ip — a p sm3ip. (11) 

This form describes both electron and hole states. For 
conduction-band electrons, there are the Rashba (c p ) and 
the Dresselhaus (b p ) terms, while the a p -term exists be- 
cause of the p-cubic Dresselhaus contribution. The u p - 
term can be attributed to higher-order invariants allowed 
by symmetry For holes in the ground-state subband, 
the dp- and &p-terms exist due to the structural inversion 
asymmetry. The term containing a p cx p 3, is the one con- 
sidered in the theory of the spin Hall effect for holes, this 
term is derived [10] from the isotropic Luttinger Hamilto- 
nian. The anisotropy of the Luttinger Hamiltonian, de- 
scribed by the parameter \i = (72 — 73)7(72 +73), where 



In application to conduction-band electrons, when the 
Dresselhaus model implies a p = Ap 3 , b p — fip — Ap 3 , and 
A ~ (3(d/2irh) 2 (for a deep square well of width d), this 
means that abruptly jumps to the universal value 
3e/8-rrh if the electron density n 2 D = Ppl'ZiT'h 2 increases 
and exceeds ir/d 2 . A similar behavior, though without a 
qualitative explanation, has been found in Rcf. 26. For 
holes, a p = —up 3 , b p = —^inp 3 , and \a\ > \b\ since \^i\ < 1. 
This means that Y,h of 2D holes in [001]-grown wells is 
insensitive to the anisotropy of the Luttinger Hamilto- 
md stays at the universal value for the case of clean 
hole systems. 

If the Cp-term is added into consideration, the anal- 
ysis leads to the phase diagram shown in Fig. 1. The 
spin-Hall conductivity is equal to 3e/87r?i at —1 < b/a < 
1 — \c/a\, in the region inside the bold triangle in Fig. 
1. There are 5 regions, and several topological transi- 
tions can take place as the parameters are varied. To 
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FIG. 3: Phase diagrams for the SOI of Eq. (11). The Berry 
phases $ for each region are indicated. The spin-Hall conduc- 
tivity is ±3e/87r7l in the regions with $ = ±37r. 

demonstrate a possibility of their experimental observa- 
tion, one should put c p = ap for electrons and c p = a^p 
for holes. The Rashba coefficient a is determined by 
structural asymmetry, while ah — 2Sa v /(hd) 2 [24], where 
Sa v = —0.035 eV nm 3 for GaAs. The results of cal- 
culations are shown in Fig. 2. For electrons, £# is 
plotted as a function of the dimensionless parameter 
2d 2 n2D/ft = (pFd/irh) 2 in the range pf < V^^U/d, when 
only the lowest electron subband in the deep square well 
is populated. If Rashba coupling is nonzero, this depen- 
dence has two jumps and the region of universal behavior 
is shifted towards higher densities. If \a/(3\ exceeds 1, 
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